arXiv:1503.00178vl [cond-mat.mes-hall] 28 Feb 2015 


Environment Assisted Quantum Transport in Organic Molecules 

Gabor Vattay,*" and Istvan Csabai,® 


One of the new discoveries in quantum biology is the role of Environment Assisted Quantum Transport (ENAQT) in excitonic 
transport processes. In disordered quantum systems transport is most efficient when the environment just destroys quantum 
interferences responsible for localization, but the coupling does not drive the system to fully classical thermal diffusion yet. This 
poised realm between the pure quantum and the semi-classical domains has not been considered in other biological transport 
processes, such as charge transport through organic molecules. Binding in receptor-ligand complexes is assumed to be static as 
electrons are assumed to be not able to cross the ligand molecule. We show that ENAQT makes cross ligand transport possible 
and efficient between certain atoms opening the way for the reorganization of the charge distribution on the receptor when the 
ligand molecule docks. This new effect can potentially change our understanding how receptors work. We demonstrate room 
temperature ENAQT on the caffeine molecule. 


1 Introduction 

There is overwhelming evidence that quantum coherence 
pWs an essential role in exciton transport in photosynthe- 
sisliHU Qjjg Qf fjjg fundamental quantum effects there is the 
Environmentally Assisted Quantum Transport®^ (ENAQT). 
It applies to partially coherent quantum transport in disordered 
systems. At low temperatures transport is dominated by quan¬ 
tum walk of excitons over the excitable sites forming a net¬ 
work. While a classical walker diffuses away from its initial 
position like ~ in time via taking random turns, a quantum 
walker takes a quantum superposition of amplitudes of alter¬ 
native paths. In a strongly disordered system the interference 
is destructive and the walker becomes stuck or ’localized’® 
At medium temperatures coupling to the environment partially 
destroys quantum interference and the walker is free to move 
and diffuse. Then at very high temperatures decoherence be¬ 
comes very distractive and the exciton gets frozen due to the 
Zeno effect®. As a result, transport is most efficient at medium 
temperatures or at medium level of decoherence and much less 
efficient at low or high temperatures. Transport efficiency is 
the highest and mean transport time is the shortest at medium 
temperatures. In photosynthetic systems parameters are such 
that this optimum is near room temperature (290K). 

The conditions for ENAQT look quite generic and one can 
suspect that it can occur in a wide range quantum processes 
in biology. Yet, the presence of ENAQT has only been estab¬ 
lished in the exciton transport of light harvesting systems. In 
this paper we demonstrate that - indeed - ENAQT can play 
a major role in electron transport through organic molecules 
at ambient temperatures. Accordingly in certain cases our un¬ 
derstanding of charge transport in biology needs revision. 
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Fig. 1 Caffeine electronic energy level spacings as a 

function of the energy £„ are indicated by red circles and connected 
with green line on a semi-logarithmic plot. Horizontal lines 
represent the mean level spacing at A = 1 AeV and the room 
temperature at hgT = 0.025eV. Energies are in the 

-35.3907eV-h57.7272eV range, the HOMO is at 

Ehomo = -11 •907eV' and the LUMO is at Eujmo = -8.6776eT. 
Calculation of A = 66 electronic levels has been carried out with the 
Extended Hiickel method. 


Charge transport in biological molecules has been studied 
extensively with a wide range of quantum chemical meth¬ 
ods®. These methods cover fully quantum, mixed quantum- 
classical, semi-classical, and fully classical approaches. How¬ 
ever, none of these cover the parameter range of the valid¬ 
ity of ENAQT. In general, electronic levels are broadened due 
to the coupling to the environment. The width of the levels 
r is proportional to the level of decoherence, which is then 
ultimately determined by the temperature of the environment 
r ^ ksT. Quantum description is relevant, when the broad¬ 
ening is much smaller than the spacing between the consec¬ 
utive energy levels E„+i —E„^Y and the mean level spac- 









ing A = {En+\ —E„) is much greater than the temperature 
A ^ hgT. The semiclassical and classical or mixed descrip¬ 
tion is relevant in the opposite limit, when A^ksT and many 
energy levels are involved in the process. The most efficient 
ENAQT sets in when the environmental temperature is com¬ 
parable with the main level spacing Ar^ hgT or at least the 
level broadening is comparableCSl with the spacings of some 
levels En+i — En ~ T ^ kgT. In small and medium sized 
(less than 500 Da) organic molecules the mean level spac¬ 
ing of electronic levels is on the electron volt scale, wich is 
about two orders of magnitude larger than the room temper¬ 
ature (about 25 milli-electron volts). So, at first sight we 
would expect that EQNAT is relevant only at temperatures 
T — leV/ks « 12000/r, which is twice of the temperature of 
the surface of Sun. However, there is a second possibility: 
In Eig. we show the nearest neighbor spacings of the en¬ 
ergy levels of caffeine. The mean level spacing is A = 1 AeV, 
which is much larger than the energy scale of the decoherence 
kgT = 0.025ey. However, there are several level spacings in 
the spectrum, which are close to 0.025ey indicating that while 
global ENAQT extending for the entire molecule is not possi¬ 
ble, there might be local pockets on the molecule, where it 
plays an important role. Next we continue with showing just 
that. 

2 Transport in ligand-receptor complexes 

One of the most significant building blocks of biochemical 
processes in the cell is the docking of signaling molecules 
in receptors. A typical situation is shown in Fig. where 
adenosine is docked in the A2A adenosine receptor. Besides 
hydrophobic interactions binding the ligand molecule to the 
receptor electrostatically, there are also hydrogen bonds con¬ 
necting the two systems. The hydrogen bridge makes possi¬ 
ble the exchange of electrons between the two systems. The 
receptor consists of amino acids and these can have very dif¬ 
ferent charging situations. Arginine, histidine and lysine have 
positive side chains while glutamic acid and asparatic acid has 
negative side chains. The rest of the amino acids is neutral 
or hydrophobic. An electron from a negatively charged part 
of the receptor protein can jump via the hydrogen bridge to 
the neighboring atom of the ligand molecule. If the ligand 
molecule would be a good conductor the electron could walk 
trough the molecule and reach a positively charged part of 
the receptor trough and other hydrogen bridge. The ligand 
molecule would act as a molecular wire, a lightning rod. As 
a result, the charge distribution on the receptor protein would 
change suddenly upon the docking of the ligand molecule and 
could spark sudden motion of the parts of the protein, due 
to the change of the equilibrium of the electrostatic forces. 
Our understanding before ENAQT ruled out such transport 
of electrons through the ligand molecule. Quantum calcula- 



Fig. 2 Schematic picture of Adenosine in complex with its A2A 
receptor. Legend: black dashed lines - hydrogen bonds; green solid 
lines - hydrophobic interactions; green dashed lines - Pi-Pi, Pi-cation 
interactions. This figure has been downloaded from RCSB Protein 
Data Bank and has been created by the Pose View softwareElE^. 

tions show that, unless tunneling plays a role in some special 
cases, ligand molecules act as insulators. Next, we show that 
ENAQT changes this picture and provides a mechanism for 
the easy and effective transport of electrons between certain 
atoms in the ligand molecule. 

3 Transport model 

We can develop a model for the description of the electron 
transfer process of the previous section. For the description of 
the molecule only the n electrons in the N atomic orbitals are 
considered. The molecular orbitals \j/j are linear combinations 
of the atomic orbitals 

N 

\j/j ^Y,Cjr(Pr and 7 ,r = 1,...,A, (1) 

(=1 

where (p,- are the valence atomic orbitals 2S,2Px,2Py and 2P^ 
for carbon atoms and hetero atoms and IS for hydrogen atoms. 
The molecular orbitals can be determined from the overlap 
matrix Srs = {(Pr \ (Ps) and Hamiltonian //„ = {(Pr \ Hg/f \ (Ps) 
matrices via the generalized eigenequation 

HC=£SC, (2) 

where H, S and C are square matrices containing the elements 
of Hrs, Srs and Cjr respectively. The raw molecular orbitals 
can be orthonormalized via the Lowdin transformation^!. The 








coefficients in the Lowdin basis can be introduced via C = 
S-i/2c' and the transformed Hamiltonian H' = 
satisfies the eigenequation 

H'C'=£C', (3) 


The time evolution is then a Slater determinant again, com¬ 
posed of for the occupied orbitals and 'T(f). 

The reduced density matrix of this wave function contains 
the occupied states and the unoccupied part 'T". In atomic 
orbital basis 


where the eigenenergies E remain unchanged. Next, we will 
drop the prime signs and use this basis as default. 

We assume that the electron is coming trough the H-bridge 
or via other mechanisms and initially enters one of the atomic 
orbits of the molecule. Initially the n electrons are placed in 
pairs on the lowest molecular orbits with opposite spins and 
occupy about the half of the orbitals. The incoming elec¬ 
tron occupies one of the atomic orbitals. In a pure quan¬ 
tum description the initial wave function is a n -f 1 dimen¬ 
sional Slater determinant of the n molecular orbitals and the 
single initial electron on the atomic orbital with wave func¬ 
tion ‘T(f = 0). This wave function is localized on the atomic 
orbital and it is proportional with the initial atomic orbital 
'T(O) = a(pi, where a is a normalization constant. The ini¬ 
tial atomic orbital can be expanded in terms of all molecular 
orbitals <P/ = L/ where is the inverse of the 

coefficient matrix. The function (pj is not orthogonal to the 
n occupied molecular orbitals and only the expansion coef¬ 
ficients corresponding to unoccupied molecular orbitals will 
contribute to the norm of the n+l electron wave function. We 
can split‘T into two orthogonal parts 

XJ/ _ IJ/O 


p=|‘P“)(‘T“|+Llrr)(rtl- 

J 

This is a very useful expression from the point of view of gen¬ 
eralization. The first term represents the evolution of the in¬ 
coming electron and the second term represents the rest of the 
electrons frozen in the Fermi sea. The incoming electron even¬ 
tually leaves the system. Just like in case of light harvesting 
systems this can be modeled by adding an imaginary sink term 
to the Hamiltonian 

H' = H-iK-1 (Pf){(Pf I, 

where K is the rate of the leaking out of the electron via a 
H-bond coupled to one of the final atomic orbital (pF- Once 
the excess electron is leaked out the reduced density matrix 
reduces to the frozen Fermi sea contribution 

P =L I V^;)(V^; I • 

j 

In our non-interacting electron approximation the leaking 
out of the electron is fully described by the non-unitary 
Schrodinger equation 


where is spanned over the occupied and ‘P" over the unoc¬ 
cupied orbitals of the molecule wit n electrons. The normal¬ 
ization condition of the Slater determinant then involves the 
unoccupied part only | ‘P" p= 1, which yields the following 
normalization condition 

L K% = ^- (4) 

y'Gunoccupied 
Consequently the norm 

\^'\^=a^ = 


^;Gunoccupied P 


ij 


is larger than unity. This reflects the fact that in this descrip¬ 
tion the incoming electron creates both electronic states in the 
unoccupied sector and hole states in the occupied sector. 

The time evolution in this model is very simple. The molec¬ 
ular orbitals are eigenfunctions of the Hamiltonian and remain 
unchanged except the stationary phase factors The 

time evolution of the incoming electron is governed by the 
Schrodinger equation 


d 

>p = H'P. 

dt 


(5) 


Pj^‘P = H'‘P, (6) 

at 

involving now the anti-Hermitian part —iPiK \ (Pf){(Pf \ de¬ 
scribing leaking. The density matrix corresponding to the 
wave function can be groupped into four sectors 

I vp)(vp 1 = 1 vp“)(‘P‘‘ I + I I + I q/")(vp« I + I vp«)(‘P“ I, 

and only the (m, u) sector has physical meaning. The evolution 
of density matrix is described by the von Neumann equation 

[H',p] =[H,p]-i{Hi,p}, (7) 

where Hi = hK \ (Pf){(Pf \ is. the leaking term and {,} denotes 
the anti-commutator. 

We can now take into account the effect of coupling of 
the electrons to the environment. This can come from many 
factors including the phonon vibrations of the molecule and 
also very crude interactions with water molecules and other 
sources of fluctuating electrostatic forces in the complicated 
biological environment of a cell. Since all the sources of 
phase breaking and dissipation cannot be accounted for we 
can treat the system statistically and use the phenomenological 





approach and add the Lindblad operator to the von Neumann 
equation 

ih^^[H\p]+J^ip). (8) 

The most general Lindblad operatoiQ^in our case can be writ¬ 
ten in terms of the projection operators of atomic orbitals 
Ar =1 (Pr}(<Pr | aS 

(2Ai.pAs AgArP pAgAr), 

rs 

where L„ is a positive definite covariance matrix of the noise 
at different atomic sites. We can get L„ from detailed models. 
The crudest approximation is when we assume a completely 
uncorrelated external noise and set where T is the 

strength of the decoherence. Its detailed form is model de¬ 
pendent, but it is in the order of thermal energy and without a 
detailed knowledge of the system can be set to r = ksT^. 

Now the steps of the full calculations can be summarized 
like this. First we identify the initial I and final F atomic 
orbitals, where the electron enters and exits the system. We 
initialize the density matrix with the initial wave function p = 
a} I v//)(v// |. The trace of this density matrix is Trp — 
and contains both occupied and unoccupied states. We evolve 
this density matrix according to (j^. The physically relevant 
reduced density matrix is given by p(“>“) the projection of the 
density matrix for the unoccupied sector and we use this for 
the calculation of physical quantities. 

4 Transit time between atoms 

The next step is to calculate physical quantities characterizing 
the electron transport between various atoms. In our model 
we have to characterize the event of a single electron passing 
trough the system. Once the electron enters the system it can 
leave the system only trough the exit. Unless some symmetry 
consideration forbids it the electron will leave the system with 
probability 1, so we have to find another characteristics. The 
next quantity is the average time T an electron needs to get 
from the entrance to the exit. This depends on two factors in 
our model; the details of the molecule and the value of K. It 
is hard to set the value of K without detailed models or mea¬ 
surements of the system, therefore we have to find its value 
from other considerations. Since the typical energy range in 
the system is the mean level spacing we can expect that the 
rate of leaking is in the order of the time (C ~ A//j. Once the 
scale of the escape rate is set we determine the average time 
needed from an initial atom to reach the destination Xpi and 
also the time an electron would need to leave the system if it 
would be placed immediately to the final site Tff- The dif¬ 
ference Xpi — '^FF is a better characteristics of the transport as 
the pure pure pumping time is transformed out and mainly the 
inter-atomic travel time enters. 


The average transit time can be calculated from the evolu¬ 
tion of the reduced density matrix. The outflowing probability 
current at the exit site F at time t is given by dP{t) = Jp {t)dt = 
2dtK{(pp I p(f) I (pp) and the average transit time is the aver¬ 
age 


POO 

T = 2Kh t {(pp \ p{t) \ (pp)dt. 

Jo 


(9) 


This integral can be calculated analytically using the solution 
of the evolution equation. Eqation (j^ can be cast into the form 


^ = ( 10 ) 

" pq 


We can re-index this equation and introduce a new index re¬ 
placing the pairs of indexes J = {r,s) and J' = {p,q). With the 
new indexing ( [T0| reads 

t = no 

The solution of this is 

p(f)=exp(^-^^p(0). 


Since the initial density matrix has only one nonzero element 
J[ = (/,/) the integral in (j^ yields the matrix element 


T = —2Kh^ 


R 


- JfPi 


( 12 ) 


where Jp = {F,F). Next, we show the results of the calcula¬ 
tion for a molecule with biological relevance. 


5 Caffeine 

In order to show ENAQT in organic molecules we should pick 
a biologically relevant example which is also computationally 
feasible. Our choice is caffeine as it is a sufficiently small 
molecule with 24 atoms (see in Table [^1. Using the Extended 
Hiickel method^^ implemented in the YAeHMOpESI software 
package. There are A = 66 valence atomic orbitals. The com¬ 
putational task involves the repeated calculation of the inverse 
of R, which is an (N x N) x {N x N) = x matrix, or 
4356 X 4356 in case of caffeine. The rapid growth of the com¬ 
putation time with is the strongest numerical limitation in 
the problem. In the calculation the escape rate K = h/leV has 
been used. There are 2145 possible pairs of initial and final 
atomic orbitals. In Eigurej^we collected the most interesting 
results. These are the transit times of electrons ending on the 
2S orbital of oxygen no. 6 in Table [T] The blue curve shows 
the transit time between this orbital and itself This is the mini¬ 
mum time Tpp needed for an electron to leak out from this site. 
This value is constant for low dephasing and starts increasing 







sible. 



Fig. 3 Transit times between atomic orbitals in the caffeine 
molecule at escape rate K" = fi/leV. Dephasing rate F in units of eV 
is shown on the horizontal axis. Curves show transition times 
between the N = 66 atomic orbitals and the 2S orbital of oxygen 
atom no. 6 in Table[T] The blue curve corresponds to the transit time 
between 25 orbital of oxygen no. 6 and itself. The green curves 
show the transit time between the 25 orbital of oxygen no. 6 and the 
orbits 2Px,2Py and 2P^ of oxygen no. 6. The red curve shows the 
transition time between the 25 orbital of oxygen no. 11 and 25 of 
oxygen no. 6. The black curves show the rest of the transit time 
curves. 


in the classical limit of large dephasing. The green curves 
show the leaking time from the same atom but the electron is 
initially placed on the three 2P orbitals of the atom. In this 
case we can observe a shallow minimum in the T = 1 — IQeV 
range. The black curves represent the transitions between all 
the atomic orbitals and the 2S orbital of oxygen no. 6, except 
2S of oxygen no. 11, which is shown in red. The black and the 
green curves show the sign of ENAQT and are large both for 
small and large T with a minimum in the middle. The mini¬ 
mum corresponds again to about 1 — IQeV in accordance with 
the expectation based on the mean level spacing lAeV, which 
is 2-3 orders of magnitude larger than room temperature. The 
red curve also has a minimum in this range. However, a sec¬ 
ond shallow minimum is also developed around O.OleV and at 
room temperature 0.025ey the transition time is still close to 
this minimum. This second optimum of ENAQT signals fast 
transport between the 2S orbits of the two oxygen atoms. 

The transition time between the two oxygen atoms is just 
10 times as large as the minimum leaking time set by the 
scale of K. This should be contrasted with the quantum limit, 
where the transition times are 1000 to 100000 times larger 
than the minimum leaking time. A pure quantum calculation 
would find the caffeine molecule strongly insulating between 
the atoms and electron transport would be practically impos- 


Table 1 Index and position of atoms in the caffeine molecule 


Index 

Atom 

X 

y 

Z 

1 

N 

1.047000 

-0.000000 

-1.312000 

2 

C 

-0.208000 

-0.000000 

-1.790000 

3 

C 

2.176000 

0.000000 

-2.246000 

4 

c 

1.285000 

-0.001000 

0.016000 

5 

N 

-1.276000 

-0.000000 

-0.971000 

6 

0 

-0.384000 

0.000000 

-2.993000 

7 

c 

-2.629000 

-0.000000 

-1.533000 

8 

c 

-1.098000 

-0.000000 

0.402000 

9 

c 

0.193000 

0.005000 

0.911000 

10 

N 

-1.934000 

-0.000000 

1.444000 

11 

0 

2.428000 

-0.000000 

0.437000 

12 

N 

0.068000 

-0.000000 

2.286000 

13 

C 

-1.251000 

-0.000000 

2.560000 

14 

c 

1.161000 

-0.000000 

3.261000 

15 

H 

1.800000 

0.001000 

-3.269000 

16 

H 

2.783000 

0.890000 

-2.082000 

17 

H 

2.783000 

-0.889000 

-2.083000 

18 

H 

-2.570000 

-0.000000 

-2.622000 

19 

H 

-3.162000 

-0.890000 

-1.198000 

20 

H 

-3.162000 

0.889000 

-1.198000 

21 

H 

-1.679000 

0.000000 

3.552000 

22 

H 

1.432000 

-1.028000 

3.503000 

23 

H 

2.024000 

0.513000 

2.839000 

24 

H 

0.839000 

0.513000 

4.167000 


6 Conclusions 

We investigated the possibility of ENAQT in organic 
molecules at room temperature. We worked out a frame¬ 
work which makes it possible to treat this problem in first ap¬ 
proximation. We showed that global ENAQT is not present 
in these molecules in general as the optimal temperature is 
several orders of magnitude higher than the room tempera¬ 
ture. However, room temperature ENAQT can be observed 
in the electron transport between certain atoms in the caffeine 
molecule. This indicates that in ligand-receptor systems the 
ligand molecule can act as a molecular wire connecting differ¬ 
ently charged areas of the receptor protein and can contribute 
to their response for the docking. This can have a signifi¬ 
cant importance in understanding signaling and drug action 
in cells. 

The authors thank D. Salahub and S. Kauffman for help, 
encouragement and numerous discussions over the last three 
years. 
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